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Abstract 



The nonlinear magnetization is considered within the Ginzburg-Landau 
theoretical framework, in the Meissner regime. A calculational method in the 
case of conventional superconductors (one order parameter) is developed and 
this method is extended for the case of two order parameters (s+d mixing). 
It is confirmed that corrections to the penetration depth, in the mean field 
analysis, are of the order of Hq^ where Hq is the applied field. We analyze 
carefully the possible solutions which lead to different scenarios in the physics 
of the symmetry of the order parameter. The anisotropy in the penetration 
depth is calculated and the temperature dependence of the magnetization is 
extracted. We discuss the relevant experimental results in the light of these 
calculations. 

PACS numbers: 74.20.De 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The symmetry of the order parameter in high-T^ superconductors is an important and 
vital question. Knowledge of the symmetry of the superconducting gap function adds in- 
formation regarding the microscopic mechanisms and reveals novel phenomena as well. An 
important contribution to this effort comes from one particular experimental technique, the 
measurement of the penetration depth or the magnetization when magnetic field is applied. 
The dependence of the penetration depth on the applied field (for small values of the field) 
and on the temperature is able to provide some conclusions on the order parameter symme- 
try. The quantity which is most easily measured is the deviation of the penetration depth 
from its zero-temperture, zero-field value: A\ab{H,T) = \ab{H,T) — Aa6(0,0). The indices 
indicate an in-plane penetration depth. The measurements that have been performed on 
different materials do not give a completely clear picture. Measurements on Y Ba2Cu'iOy 
crystals by Sridhar et al. [|l| result in a quadratic field-dependent penetration depth below 
Hci and linear dependence above Hd. On the other hand Bi2Sr2CaCu20y (BSCCO) shows 
a quadratic behavior in a region of values of magnetic field much smaller than Hd and 
linear dependence in the remaining region 0. The temperature dependence shows a very 
interesting behavior as well. Hardy et al. observed linear T-dependence at low temperatures 
in YBCO and Bonn et al. observed a crossover from linear-T to upon Zn doping 
0]. Also a behavior has been observed on BSCCO by Maeda et al. f^. The linear 
field-dependence is predicted by the d-wave scenario, according to Xu et al. ^j. 

At this moment the question of the microscopic mechanism is still lacking a clear con- 
sensus, though most work now concentrates on the magnetic mechanism (antiferromagnetic 
fluctuations in Cu02 planes) as the first candidate for the d-wave component. Experimental 
results support an order parameter of predominantly d-wave symmetry, as reviewed in Refs. 
and 1^. However, there are possible indications of an admixture of s-wave in overdoped 
and electron-doped materials, as reviewed in Ref. 0. This latter possibility is consistent 
with the magnetic hypothesis. 
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A signature of s-d mixing is the spontaneous appearance of orthorhombic-like anisotropies 
in the superconducting state. In the penetration depth this would mean A^. 7^ Aj^. In a system 
such as YBCO with an orthorhombic crystal structure, it is necessary to somehow disentangle 
this spontaneous anisotropy from the crystalline anisotropy. This paper attempts to do this 
for one particular measurable quantity, the temperature- and field-dependent penetration 
depth. 

Experiments in YBCO have indeed shown a large difference in the penetration depth 
along the a axis (Aa) and b axis (A;,), which has been attributed to the orthorhombic distor- 



tion |10,11|. This distortion can be taken into account in the construction of the free energy. 



In Ref. the anisotropy at zero temperature has been calculated in a microscopic model 
where the pairing is of d-wave type and using a single tight-binding band with different hop- 
ping parameters in the a and b directions. The point of view is that the chains accomodate 
a significant part of the condensate and this has to be taken into account. Also in Ref. 
p!3| it was found that the chains cannot become superconducting by proximity effect. The 
microscopic theory of the non-linear Meissner effect within a purely one-component d-wave 
scenario has been formulated by Xu et al. Their emphasis was on low temperature 
behavior, while our Ginzburg-Landau approach is more suited to temperatures neat T^. 

In this paper, the framework of Ginzburg-Landau (G.L.) theory is used and this mean 
field theory allows us to do calculations taking into consideration the possible symmetry 
of the order parameter and constructing the appropriate free energy functional. In this 
manner the specific microscopic mechanism is avoided and the thermodynamic properties 
can be extracted. The attention is concentrated in the Meissner phase, where the presence 
of vortices can be neglected and the field can be considered weak. The solution of the G.L. 
equations beyond the London approximation offers the dependence of the penetration depth 
and the magnetization on low applied fields and also the temperature. In mean field theory 
then the dependence on the field is quadratic. The corrections of the order parameter can 
be calculated analytically as well. 

The plan of the paper is as follows. In Sec. II, the G.L. equations are solved in the 



case of conventional superconductors in the Meissner phase, so the calculational method is 
illustrated. We believe it is really worth presenting these calculations, since we were unable 
to find any detailed calculation in the literature. In Sec. Ill, the method is extended to the 
case of the mixed symmetry in the order parameter ( s + d). The possible solutions are 
considered and we discuss the implication of each one. In general the situation where the 
d-wave is dominant is adopted and the role of the s-wave component is illustrated. This 
s-symmetry part appears partly because of the perturbation due to the field but a nonzero 
value even in the absence of external field can be present. The anisotropy in the penetration 
length along the a and b axes can be calculated easily since the orthorhombic distortion is 
considered in the development of the free energy functional. In Sec. IV we conclude and 
discuss the results in connection to relevant experimental data. Also in the Appendix the full 
set of equations are given and the implication of the different field directions are considered. 



II. ONE ORDER PARAMETER - CONVENTIONAL SUPERCONDUCTORS 



For the case of the conventional (s-wave) superconductors, the G.L. free energy is written: 



a|$(x)p + /3|^(x)|^ + 7|n*|2 + ^ 

ovr 



(1) 



where 11 = —ih'V—'^A and i? = V x A. The above free-energy is minimized with respect 
to \I'* and A. ^ is written as ^ = |\I'| e**^, where is the phase of the order parameter. 
Then the equations to be solved are : 

2e 
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+ a^ + 2/3|^|^^ = 



c ^ ^ ^ ^ 2e ^ 

— V X V X A = 47e|^^P(W0 A) 

An c 



(2) 
(3) 



In the London approximation is considered spatially constant, unaltered by the ap- 



plied field. From Eq. (]^) above the penetration depth is found to be A = \J {—c^jS/lGn'yae). 
In the calculations, the applied magnetic field Hq is taken parallel to the z-axis, the super- 
conductor occupies the half space x > and the most convenient gauge in which to work 



is the one with Ay = Az = 0. Then the phase and the vector potential A^ have x,y- de- 
pendence only. Since the field and the supercurrent are expected to have only x-dependence 
(translational invariance in the y direction) (p and A^ can be written as : (p = y g{x) and 
Ax = y a{x). is expected to have x-dependence as well. After these substitutions (also 
for convenience |\1'| = /) the above equations become : 

7(-^'0 + ^V/) + 2/?/^ + «/ = O (4) 

= ^^'^^ 

Eq. @/® above is the real/imaginary part of Eq. (H) respectively. The boundary 
conditions are: a(x = 0) = Hq, a{x = oo) = 0, g'{x = 0) = {hc/2e)HQ, g'{x = oo) = 0, 



f{x = oo) = /o = \J —a/213 and f'{x = 0) = IQ- It is obvious then that g oc Hq and the 
second term in the l.h.s. of Eq. (^) above serves as a term containing the small parameter 
of the problem. Then perturbation theory can be used to obtain the solution in small fields 
(small compared to the thermodynamic field He). We consider solutions of the form: 

f = f, + Hlh + H',h + ... (7) 

In this case, if the above ansatz substitutes f in Eq. (^ and terms with same power in 
Hq are gathered {g ^ g / Hq) we get: 

-7h'^ + 2Pf^ + afQ = (8) 
-7^'0 + + «/i = -ih'g'fQ (9) 



he d g 2eh 



{fl + 2Hlhh) ~g (10) 



Svre dx'^ 

The last equation comes from the combination of Eqs. (^ and (P). The above expansion 
is correct, since the equations corresponding to greater powers of Hq are such that conver- 
gence is guaranteed. Therefore the above method works and leads to correct results. Then 
the largest correction due to the field will have an Hq dependence. The solution for /o is 



the unperturbed part y— a/2/9. The solution for /i then can be found easily and the form 
of the order parameter becomes (He = ^^^J 



m = ^J-a/2(3{l + A 



2 

^ exp(-x/V20 + exp(-2x/A) 



A 

An interesting immediate observation is that besides the term with the exponential 
dependence on the coherence length there is another one which decays with the characteristic 
length of A/2. The above values of ^ and A are the zero- field ones. In the case of a 
superconductor with ^ = \/2^/2 the divergent part in the denominator is cancelled by 
the numerator. Also, the nonlinear correction /i is negative, independent of the type of the 
superconductor, something which is expected since the magnetic field acts as a pair-breaking 
mechanism. 

The next step is to substitute the above solution in the equation of g and solve it by the 
same method. We define the effective penetration depth as follows: 

1 f°° 

^eff = -jT / a{x)dx (12) 
By performing the integration the effective penetration depth becomes : 

The penetration field is an increasing function of the field and has a quadratic dependence 
on it. It is somewhat remarkable that this expression does not appear anywhere in the 
published literature, as far as we are aware. It holds for every value of n = X{H = 0)/^{H = 
0), and is therefore a generalization of a formula given in Ref. 0: 

<5A/A = 3/4 X {Ho/H.Y, (14) 

which holds in the limit /t — oo. The finite k correction is due to the part of the order 
parameter that decays with the characteristic length of A/2. 
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III. S+D ORDER PARAMETER SYMMETRY 



Recent experimental evidence and also theoretical calculations |15] suggest that a mixed- 
symmetry order parameter is a possible candidate to explain several features observed in 
experiments on high-Tc materials in the absence of a field. However, even if there is no 
s-wave component in zero field, an s-wave component is always formed in the vicinity of 
vortex cores or due to induced currents or due to surface or impurity effects or, finally, 
as a result of the orthorhombic distortion in these materials. Within the Ginzburg-Landau 
theoretical framework, the free-energy functional that takes into account both the symmetry 
of the material and the order parameter (but with fourth-order derivatives neglected) takes 



the form ||16| : 



F = j (fx [aa\^d? + md\' + d? + as\^s? + Al^/ + Is^s? + /Jsl^sH^dr (15) 

+ /34(*s*'^/ + ^/'^s') + 7.[(n,^.)*(n,^,) - (n,^,)*(n,^,) + c.c] + ^| . (i6) 

The orthorhombic distortion has been taken into account, through the mixed-gradient 
term (and a perturbative bilinear on the two components term). These terms may serve as 
a "source" for the s-component (it's one of the two "scenaria" that come from the possible 
solutions). The mixed-gradient term distinguishes a-axis from b-axis as well (c-axis is along 
the z-direction) . The above functional can be used to derive the anisotropy in the penetration 
depth that is observed. The same G.L. functional has been extensively studied recently 
for the case of vortices and actually has been derived in the weak-coupling limit for both 
continuous and lattice Hamiltonians [|r^-pO|. The geometry will be the same (Meissner- 



geometry) and the difference is that now penetration along both x and y directions have 
to be considered separately. Another difference comes from the two-dimensional nature of 
the free energy functional. The orientation of the applied field produces somewhat different 
results so one has to consider both the "in-plane" case and the case where the field is along 
the z-direction. 
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A. London Approximation 

After minimization Eq. [T^ with respect to "^d* and A the Euler-Lagrange equations 
to be solved are: 

i-fdU^ + ad)^d + 7.(n/ - n,')^. + 2(32\^d\^^d + /?3|^.r^. + 2/?4^.'^/ = (17) 

(7.n2 + as)-^s + 7.(n/ - n,')^, + 2A|^,|2^, + Ps\^d\^^d + 2/34^/*/ = (is) 
X V X 1 = ^ {7d*/n^d + 7.^s*n^d + 7.(^/(^n, - ^n,)^, + 

^/(yn,-xn,.)^s + c.c.} (19) 

From the form of the above two first equations it's easy to verify that the only possibility 
for the relative phase of the two components is or vr, therefore only d±s states are 
the starting point of our analysis. Physically d+s and d-s are equivalent and the system 
spontaneously chooses one of these states. The two axes x and y are identified with the 
crystallographic a and b in order to make connections with the experiments. The applied 
field is along the ^-direction and its spatial variation is along the x-direction. The boundary 
conditions are: 

|\E's(x = oo)| = fso = const. l^^dix = oo)| = fdo = const. A{x = oo) = 

[-fd^^d + iviyTiy - xU^)^s] ■x = o [-fsii^s + iviyT^y - £n^)^d] ■x = o 

IvE'sol = fso and |\l'do| = fdo are the bulk values of the two order parameters, without 
the applied field. If the London approximation is made (|\l/d|, l^^sl spatially constant) the 
penetration depth along the two axes can be found to be: 

\ ^ r_r 11/2 C2^^ 

^ Unei^dl^dl' + lsl^sl'+lvl^sW^dW ^ ' 

r2 1 
\ ^ r_r 11/2 

' U7re(7d|vl/rf|2 + 7,|vI/,|2-^^|vI/,||vI;,|)J ^ > 

The above formulas give a first estimate how it is possible to obtain the difference in the 



penetration depth from the G.L. theory consistent with the experimental measurements [|10 
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The ratio \x/\ is 1 — e where e is the quantity [^°[ + 0{{^^ |§4)^)- -^^S- HI ^^^^^ ratio 



which measures the anisotropy has been plotted as function of temperature and, as we see 
close to Tg where the s-component grows as we move from the higher temperatures towards 
the lower ones, the anisotropy becomes stronger. Since the existence of the s-component 
enhances the anisotropy, close to Tc the anisotropy vanishes. 

What this calculation demonstrates is the fact that the anisotropy can shhow a strong 
signature of s-d mixing if the ratio of the two components is temperature-dependent. The 
photoemission measurements on BSCCO-2212 ||22| have been interpreted as showing just 



such a dependence p3[. This interpretation can therefore be checked by magnetization 
measurements, if careful attention is paid to anisotropy. 



B. Beyond the London Approximation 

For the calculations beyond the London Approximation the method of Sec. II can be 
followed closely. From the nature of the equations there are two points of view in obtaining 
the solutions, with substantial difference in the physics of the problem. The first point of 
view, or the first physical "scenario" is the one starting from the assumption that > 0, so 
that the zeroth order (in the magnetic field) value of the magnitude of the \E's = fso is 0. In 
other words, the onset of the s-component of the order parameter is caused by the application 
of the magnetic field and the induced currents and can be viewed as a perturbation to a 
robust d-wave superconductor, or a small "transformation" from d to d+s, close to the 
boundary. There are three characteristic length scales in the problem, the two coherence 
lengths and that characterize the spatial change of the two distinct order parameters and 
the penetration depth A which is the characteristic length of the electromagnetic changes. 

The second physical picture is when < and a nonzero value of fso at zero current is 
considered. Then the superconductivity is truly two-channel and the existence of the s-wave 
part has to do either with the pairing mechanism or with the departure from the tetragonal 
symmetry of the lattice and not with the applied magnetic field which at most modifies the 
form of the order parameter close to the boundary. In this picture there are four different 
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characteristic length scales, the two coherence lengths ^d, the two penetration depths 

Xd and Xg (combined to one A) : 

^ = A + ^±^^tV (22) 

Xx,y Xd Xg -s/ldll Xd Xs 

1. Zero \^ s\ o-i zero applied field 

In this case there are the following consequences : (i) The onset of s-wave is due to the 
mixed gradient term and the value of the order parameter l^^l oc 7^ x Hq^. (ii) Since this 
mixed gradient term is responsible for the s-wave it is clear that when the y-dependence is 
examined alone one gets |(i + s| state, on the other hand when the x-dcpcndcncc is examined 
one gets |o? — s|. (iii) According to (i) and (ii) anisotropy cannot be observed in penetration 
depth due to term proportional to 7^ in Eq.[15] and [16]. This term finally gets proportional 
to 7^^ and the sign difference cancels out. (iv) The only way to get the anisotropy in 
penetration depth is to consider different values for the "masses" in different directions. 
Namely ^ydx and ^ydy and the same for 7^ possibly, (v) Due to the above observations the 
quantity X^/ Xy is temperature- independent. 

The two coherence lengths are given by the equations: 



^'^^ (23) 



From these expressions it can be seen that if ctg = OLg^ilT — and ad — 0(do{T — Td) 
{Td > Tg) then the actual second order phase transition is at 7^. At Tg there is only a 
crossover, not accompanied by a phase transition. The coherence length doesn't diverge 
at Tg but at a temperature much closer to Td- 

2. Nonzero 1^*5 1 at zero field 



This case can be considered naturally as a consequence of orthorhombic distortion of 
the lattice. Then the two representations s and dx2^y2 can be mixed and a bilinear term 
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in and has to be included in the free energy. The calculations are presented in the 
Appendix. The main result is that now the penetration depth acquires an anisotropy due 
to the nonzero value of "^so- The important difference from the case with nonzero \I^so is 
the different temperature dependence which actually can distinguish between the two cases. 
The distortion of the order parameter is given by: 

l^dl = l^dol + (A,, e--/v^«^ + A,, + a,, e'^/^^^) (24) 

\^s\ = l^sol + {^f (Si. e-^'/^«^ + e-'^/' + Bss e'^/^^'^) (25) 

The coefficients have been calculated in the Appendix and they depend on the G.L. 
coefficients (terms beyond the first order in 7^, have been neglected). The thermodynamic 
field He has been taken as Hn = J^l , ■ For the -u-direction the results are the same with 
the substitution 7^, —■jv 

Knowing the corrections of the order parameters, the effective penetration depth can 
be computed easily (see Appendix). The main point of the calculation is the temperature 
dependence of the anisotropy Xx/\- In Fig. ^ the temperature dependence of the ratio of 
relative corrections due to the field ^^""'^^'^ is plotted, in the limit of strong type-II supercon- 
ductivity (k >> 1). This quantity is plotted in order to avoid the explicit field-dependence. 
The relative corrections are larger in the temperature regime where the anisotropy of the 
penetration depth at zero field is larger. The effect of a small term which measures the 
orthorhombicity would be to get smooth curves in the region close to T^. 

Again we see that the appearance of the s-wave admixture gives a strong signature in 
the penetration depth, this time in the nonlinear signal. 

IV. DISCUSSION - CONCLUSIONS 

The simple picture that this paper offers may be able to give a qualitative point of view 
of the several features that are observed simultaneously in penetration depth experiments. 
The basic conclusions are : 
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(i) The anisotropy in the penetration depth may arise either as a consequence of directional- 
dependent "masses" or as a consequence of the s and d-wave mixing. As the figures show, 
mixing may resuh in strong temperature and field dependence, so the cases can be distin- 



guished. In Ref. ||T0[, considerable anisotropy was observed: Xa > 1.5A;,, but the temperature 
dependence was not measured. If the anisotropy is temperature dependent then the scenario 
with the non-zero fso has to be investigated more carefully. It is ineresting that Xa/Xb ap- 
peared to be strongly sample-dependent in these experiments. While s-d mixing is expected 
to be a strong function of doping, the effective mass ratio is not. 

(ii) The field - dependence of the penetration depth is of order Hq^ for both cases at least in 
the geometry that it is described. This is in fact a consequence of the boundary conditions 
in the studied geometry. 

(iii) The temperature dependence of the anisotropy can be derived explicitly and compared 
to experiments. Of course the G.L. equations are not valid at low temperatures, therefore 
nothing can be said for the experimentally observed crossover from T to dependence of 
the quantity 1/A(T)^ §. 

(iv) Enhancement of the coupling and the transport along the c-axis if the applied field is in 
the a — b plane. Then it's unavoidable to get screening currents along the c-axis and actually 
the penetration depth has slightly different values from the ones observed when Hq/ / z. 

(v) Finally the question of d+s versus d-s states arises (and consequently the equivalence of 
the two directions x and y. In the model that is described through the G.L. equations for a 
tetragonal crystal structure a spontaneous symmetry breaking is assumed that distinguishes 
between the two states and directions. The system chooses one of the two equivalent states. 
Measurements of the effects described here become much more complicated if domains of 
d+s and d-s form. In YBCO experiments, this requires detwinned crystals. 

In Ref. the nonlinear Meissner effect was examined in the framework of a microscopic 
model. The central mechanism that produces the nonlinearity is the different contributions 
to the current from the "co-moving" (condensate velocity parallel to Fermi velocity) and 
"counter-moving" (condensate velocity antiparallel to Fermi velocity) excitations close to 
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Fermi energy. Then in the case of pure d-wave the correction is hnear in the field. But even 
at very low temperatures this model predicts a transition to a regime where the corrections 
are of the order Hq^. In addition to that a very small admixture of another component is 
able to change the behavior of the penetration depth at very low temperatures as a function 
of temperature. Our picture is complementary to this in that it is nonly appropriate for 
relatively high temperatures. The effects predicted involve relative values of the different 
gap components which are realistic. The picture has the advantage of being rather general: 
note that effects such as impurity scattering are included in the parameters of the G.L. 
functional. 

There are several remarks on some other issues of the calculation of the penetration 
depth. The most important of these is the role of different type of fluctuations. In Ref. 

has been calculated the existence of Off Diagonal Long Range Order (ODLRO) in the 
different "states" of the superconductor. It is found that in the Meissner state this phase 
coherence is destroyed by phase fluctuations below two dimensions. Also only in strongly 
type II superconductors, the fluctuations of the fleld can be considered unimportant. In Ref. 
P5| the critical fluctuations in the penetration depth as the transition is approached from 
below is studied and it is found that both the penetration depth and the coherence length 
diverge with the same exponent which has value (0.53) very close to the mean- fleld one in 
three dimensions which was obtained in experimental work of Ref. in contrast with the 
value of ~ 1/3 obtained in Ref. p6[ on YBa2CusOQ,g5 pure bicrystals. These issues suggest 



flrst that the calculations presented here are valid in a region not very close to since in 
any case the small parameter of the problem is in fact the applied fleld in comparison to 
thermodynamic fleld which in turn is temperature dependent. Second the fluctuations in 
the case of two order parameters will have to be included in a subsequent paper because of 
the sensitivity of the second order parameter which is itself a "secondary" effect. 

In summary we have presented the mean-fleld analysis in the Meissner state of a super- 
conductor in the case of one and two-component order parameter which is relevant for the 
high-Tc materials. In the latter case the picture of the onset of the second component due to 
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different reasons have been presented and several predictions and qualitative explanations 
have been made. 
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VI. APPENDIX 

We consider the case of the applied magnetic field Hq parallel to the c-axis and varying 
along the x direction. The geometry is the same as in Part II. By following the same 
considerations as in Part II the phase can be written as (f) = y g{x) and the gauge A = 
y a{x)x. The magnitudes of the order parameters are denoted as = fs and |\E'd| = fd- 
We substitude all the above into the Eqs. ([TTD and (|l^). Then by dividing the equations 



into real and imaginary parts and also by taking into account the power of y, we obtain the 
independent equations: 

= -^a{x) 26) 
ax en 



7d^'(-^ + g'fd) + + 7.^'(0 + + W2h' + = (27) 



7.^'(-0 + g'fs) + asfs - 7.^'(^ + g'h) + P'.h'fs = O (28) 

The magnetic field (which is explicitly contained in g) is considered as the small param- 
eter of the problem, as in the case of the single order parameter and we seek solutions of 
the form: 

9 = Hog (29) 
fd = fdo + Ho' U (30) 
fs = fso + Ho' U (31) 
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An important point here is that the dependence of fd and fs on the field is quadratic 
because of the boundary conditions at x = 0. With different boundary conditions, hnear 
dependence on the field can be obtained. But this is not the case here. Finally, by gathering 
terms with the same power in the applied field, the equations for fdo, fdi, fso, fsi are: 

adfdo + 2f32fdo^ = (32) 

asf so + 2/?i/.o' = (33) 

7-i^'(-^ + g'fdo) + aaU + 7.^'(^ + ffso) + 6/?2/doVdi + P'sfso'U = (34) 

%^'(-^ + ffso) + asf si - lvh\^ + ffdo) + (3'Jm'fsi = (35) 



The possible solutions of the Eq. (|33D are /s = and fs = ^J—j^^- In the results 
presented in the text the difference in the two cases is discussed. The calculations here 
are presented under the hypothesis of a nonzero fs- From Eq. (E2, Rl) the zeroth order 



value of the two order parameters can be easily obtained. Then the equations (|3^ , |35D 
suggest the following : It is natural to expect dependence of fdi on the coherence length 
as well as the penetration depth A (exactly like the case of one order parameter) but due 
to the appearance of the terms which depend on it is also expected a small dependence 
on the coherence length C,s- The same, by symmetry, is true for the correction of the s- 



component fsi- So, the homogeneous parts of Eq. (|3^ , |35|) have solutions which depend 
on the corresponding coherence length of each component. Then the non-homogeneous 
differential equations mix the dependences as discussed above. For the exact form of the 
solutions, the boundary conditions have to be taken into account. They take the explicit 
form : fsi{x = oo) = fdi{x = oo) = 0, /^(x = 0) = f'^{x = 0) = and a{x = 0) = —Hq, 
a{x = oo) = 0. Having elaborated the crucial points, it is straightforward, after some 
lengthy algebra to obtain the full solutions-corrections to the two order parameters. The 
homogeneous parts of the equations accept solutions of the form : fdi^h = ^id exp{—x/ V^Cd) 
and fsi,h = Bis exp{—x / \/2^s) ■ If we substitute these expressions into the homogeneous 
parts of the differential equations, we obtain the expressions for the two coherence lengths, 
which are the same as in the case of zero fso up to order 7^. The solutions of the non- 
15 



homogeneous parts are : fdi,inh = exp{-2x/\) + A^d exp{-x/y/2Q and fsi,inh = 
i?2s exp{—2x/\) + B-^s exp{—x/^/2^d)- The coefficients can be evaluated by substituting the 
above expressions into the inhomogeneous parts of the differential equations and by using 
the boundary conditions. Finally we get : 

, 4^2 . -7d/do(-4/A^ + 1/2^.^) + 7./so(8/A^ - 1/2^.^) 
^ A7.(-4/A2 + l/2e/)(-4/A2 + l/2e.') 

^es'7d(l-&/e.)^)(4/A2-l/2e.^)J ^ ^ 

, _ , -7d/do(-4/A^ + 1/2^.^) + 7./.o(8/A^ - 1/2^.^) . . 

e/7.(-4/A2 + l/2e/)(-4/A2 + l/2e.^) ^''^ 

~ ^ - fe/^.)^)(-4/A2 + 1/2^.^) ^''^ 

R , /7T /. 2r A(-7j.o(-4/A^ + 1/2^/) + 7./.o(8/A^ - 1/2^/)) 
^ es7s(-4/A2 + l/2e/)(-4/A2 + 1/2^.^) 

A7.(l-(We.)^)(4/A2-l/2e/)^ ^ ' 

R , -7J.o(-4/A^ + 1/2^/) + 7Jdo(8/A^ - 1/2^/) 

e/%(-4/A2 + l/2e/)(-4/A2 + l/2e.^) ^'''^ 



5„ = 4v/2 f41) 

7desed(l-fe/e.)^)(-4/A2 + l/2e.^) ^ ' 

The effective penetration depth then can be obtained from the solution of the Eq. ( pi?]) 
with the appropriate boundary conditions that have been described, taking into account Eq. 
(p6|). The calculation is straightforward and the result is : 

The coefficients i=l,2,3 above are given by : 

Ci = Aid + -Bss 

C2 = — A2d + -j^ B2S 

E F 
C3 = Am + -j^ Bis 



where D = 'jdfdo^ + Isfso^ + Ivfsofdo, E = 2'jdfdo + Ivfso and F = 27^/^0 + Ivf, 



do- 
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FIGURES 

FIG. 1. The anisotropy \x/\ at zero field as a function of temperature. The values of the 

parameters that have been used are = 0.1 {^-^zyY^"^ ^ Is/ld = 1 and 7^ = 0.67^. 

FIG. 2. The ratio of the relative corrections due to the applied field in the regime of strong 
type-II superconductor. The values of the parameters are the same as in Fig. 1. In addition to 
that Pi=h= /?y2 
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